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ABSTRACT 

A cell membrane can be simply regarded as composite material consisting of 
lipid bilayer, membrane cytoskeleton beneath lipid bilayer, and proteins embed- 
ded in lipid bilayer and linked with membrane cytoskeleton if one only concerns 
its mechanical properties. In this Chapter, above all, the authors give a brief 
introduction to some important work on mechanical properties of lipid bilay- 
ers following Hclfrich's seminal work on spontaneous curvature energy of lipid 
bilayers. Next, the entropy of a polymer confined in a curved surface and the 
free energy of membrane cytoskeleton are obtained by scaling analysis. It is 
found that the free energy of cell membranes has the form of the in-plane strain 
energy plus Helfrich's curvature energy. The equations to describe equilibrium 
shapes and in-plane strains of cell membranes by osmotic pressures are obtained 
by taking the first order variation of the total free energy containing the elas- 
tic free energy, the surface tension energy and the term induced by osmotic 
pressure. The stability of spherical cell membrane is discussed and the critical 
pressure is found to be much larger than that of spherical lipid bilayer without 
membrane cytoskeleton. Lastly, the authors try to extend the present static me- 
chanical model of cell membranes to the cell structure dynamics by proposing 
a group of coupling equations involving tensegrity architecture of cytoskeleton, 
fluid dynamics of cytoplasm and elasticities of cell membranes. 

INTRODUCTION 

Cells are the basic elements of living organisms, such as plants and animals. 
Viewed from structure, an eukaryotic cell consists of cell membrane, cytoplasm, 
cytoskeleton, nucleus and so on. Cell membrane defines the boundary between 
the living cell and its environment, including extracellular matrix and liquid 



surroundings. Cytoplasm is a kind of viscous fluid. Cytoskeleton is a three- 
dimensional structure composed of three classes of fibers [1]: microtubules (20 
nm in diameter), elements built of polymers of the protein tubulin; microfila- 
ments (7 nm in diameter), built of the protein actin; and intermediate filaments 
(10 nm in diameter), built of one or more rod-like protein subunits. Nucleus of 
the eukaryotic cell is enclosed by membrane. 

Membranes consist of lipids, proteins and carbohydrates etc. Lipids and 
proteins are dominant components of membranes. One of the principal types of 
lipids in membranes is phospholipid. A phospholipid molecule has a polar hy- 
drophilic head group and two hydrophobic hydrocarbon tails. In physical point 
of view, a lipid molecule can be simply regarded as an amphipathic rod. When 
a quantity of lipid molecules disperse in water, they will assemble themselves 
into a bilayer in which the hydrophilic heads shield the hydrophobic tails from 
the water surroundings because of the hydrophobic forces. 

There are many simplified models for cell membranes in history [2] . Among 
them, the widely accepted one is the fluid mosaic model proposed by Singer 
and Nicolson in 1972 [3]. In their model, cell membrane is considered as a lipid 
bilayer where the lipid molecules can move freely in the membrane surface like 
fluid, while the proteins are embedded in the lipid bilayer. Some proteins are 
called integral membrane proteins because they traverse entirely in the lipid 
bilayer and play the role of information and matter communications with both 
the inside of the cell and its outer environment. The others are called peripheral 
membrane proteins because they are partially embedded in the bilayer and 
accomplish the other biological functions. Beneath the lipid membrane, the 
membrane cytoskeleton, a network of proteins, links with the proteins in lipid 
membrane. 

The first step to study the elasticity of cell membrane is to study lipid bilayer. 
Usually, the thickness of lipid bilayer is about 4 nanometers which is much less 
than the scale of the whole lipid bilayer (about several micrometers). Therefore, 
it is reasonable to describe the lipid bilayer by a mathematical surface. In 
1973, Helfrich [4] recognized that the lipid bilayer was just like a nematic liquid 
crystal at room temperature. Based on the elastic theory of liquid crystal [5], 
he proposed the curvature energy per unit area of the bilayer 



where k c is an elastic constant; and H, Co are mean curvature and spontaneous 
curvature of the membrane surface, respectively. In Gaussian curvature K 
is not written explicitly because only close bilayer is discussed in this chapter. 
We can safely ignore the thermodynamic fluctuation of the curved bilayer at 
the room temperature because of k c w 10 _19 J 3> ksT [6,7], where fee is the 
Boltzmann factor and T the room temperature. Based on Helfrich's curvature 
energy, we can express the free energy of a closed bilayer under the osmotic 
pressure Ap (the outer pressure minus the inner one) as: 



f c = (fc c /2)(2iJ + c ) 2 , 



(1) 
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where dA is the area element and V the volume enclosed by the closed bilayer. 
jU is the surface tension of the bilayer. Based on above free energy many re- 
searchers studied the shapes of bilayers [8-10]. Especially by taking the first 
order variation of above free energy and doing the complicated calculations of 
tensors, Ou-Yang et al. derived an equation to describe the equilibrium shape 
of the bilayer as [11, 12]: 



Ap - 2(xH + k c (2H + c )(2H 2 - c H - 2K) + k c X7 2 (2H) = 0. (3) 



This equation is now called the shape equation of closed membranes or gen- 
eralized Laplace equation. They also obtained that the threshold pressure for 
instability of spherical bilayer was Ap c ~ k c /R 3 , where R being the radius of 
spherical bilayer. 

Using the shape equation J2J of closed bilayers, Ou-Yang predicted that there 
was a lipid torus with the ratio of two radii being exactly v2 [13] . His prediction 
was soon confirmed by the experiments [14-16] . Otherwise, researchers obtained 
an analytical solution to Eq.© which explained the classical problem [17] — the 
biconcave discoidal shape of normal red cells [18]. 

Recently, Tu and Ou-Yang have proposed a mathematical scheme to discuss 
the elasticities and stabilities of cell membranes with membrane cytoskeleton 
and found that the membrane cytoskeleton enhances the stabilities of cell mem- 
branes [19]. But they have omitted the effect of in-plane modes of membranes 
on the stabilities. In this chapter, we will first retrospect to elasticities and 
stabilities of lipid bilayers following Helfrich's seminal work on spontaneous 
curvature energy of lipid bilayers. Next, we will fully discuss the entropic elas- 
ticity of membrane cytoskeleton, as well as the elasticities and stabilities of cell 
membranes with membrane cytoskeleton. Last, we will expatiate on how to 
construct the framework of cell structure dynamics involving tensegrity archi- 
tecture of cell cytoskeleton, fluid dynamics of cytoplasm and elasticities of cell 
membranes (with membrane cytoskeleton). 



In this section, we will recur to some main results on elasticities and stabilities 
of lipid bilayers by adopting the mathematical scheme proposed in Ref. [19-21]. 

We use a smooth and closed surface M in 3-dimensional Euclid space E 3 to 
represent a membrane. As shown in Fig. ^ we can construct a right-hand or- 
thonormal system {ei, e2, at any point r in the surface and call {r; ei, &2, e^} 
a moving frame. The differential of the frame is denoted by 



where u>i, loi and w%j — —<-<Jji — 1, 2, 3) are 1-forms. The structure equations 
of the surface are 



LIPID MEMBRANE 




(4) 



du>i = W12 Aw 2 ; 



(5) 
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doj 2 = lo 2 \ A ui ; (6) 
Wi3 = owi + buj 2 , UJ23 = buii + coj 2 ; (7) 
doJij = ujikAukj = 1,2,3). (8) 

Readers should notice that the operator "<i" is an exterior differential operator 
[19] in this chapter. The area element, mean curvature and Gaussian curvature 
are respectively expressed as [22]: 

dA = uux Aoj 2 , (9) 
H = (a + c)/2, (10) 
K = ac-b 2 . (11) 



e 3 




Figure 1: Smooth and orientable surface M. we can construct a right-hand or- 
thonormal system {ei, e 2 , e^} at any point r in the surface and call {r; ei, e2, 63} 
a moving frame. 



If M undergoes an infinitesimal deformation such that every point r of M 
has a displacement St, we obtain a new surface M' — {r'|r' = r + <5r}. 5r is 
called the variation of surface M and can be expressed as 

Sr = Sir + 5 2 r + S 3 r, (12) 
S t r = a t e t (i= 1,2,3), (13) 

where the repeated subindexes do not represent Einstein summation. Due to 
the deformation of M, ei,e 2 ,e 3 also change. We denote the change as 

Siei = Slujej, fluj = —Cliji. (14) 

Using the commutativity between Si (i — 1,2,3) and d, we obtain the fun- 
damentally variational identities of the move frame [19]: 

S1U1 — dtli — (15) 



4 



U 1 UJ 2 


— O1 Ci)m — ,Oi 1 n 


(16) 


^113 


= aQi, tt 123 = btt 1 ; 


(17) 


S 2 <JJl 




(18) 


82OJ2 


= dQ 2 - wif2 2 i2, 


(19) 


^213 


= 6fi 2 , ^223 = c0 2 ; 


(20) 


S 3 UJi 


= ^3^31 — W 2 ^321, 


(21) 


5 3 U> 2 


= ^3^32 ^ ^1^312, 


(22) 


dn 3 ■- 


= ^313^1 + ^3230;2; 


(23) 


SlLUij 


— dfllij + QukLOkj — UJik^llkj- 


(24) 



Now we can take the first order variation of functional © by considering 
above variational identities 1|15[1 - (|24[1 and Stokes theorem. It is not hard to 
obtain [19]: 

8T H = J [Ap - 2 f iH + k c (2H + c )(2H 2 - c H - 2K) 

+ k c y 2 (2H)]Q 3 dA. (25) 

Thus the Euler-Lagrange equation corresponding to the functional Th reduces 
to equation JJJ) because fl 3 is an arbitrary function. Equation Q describes the 
equilibrium shapes of lipid membranes. 

Similarly we can obtain the second order variation of functional (J2J : 

8 2 T H = J ^sV 2 ^ 3 [k c (UH 2 + 2Hc Q - 4K - c^/2) ~ n]dA + k c J (V 2 fl 3 ) 2 dA 

+ Jk c (2H + c ) S7{2Hn 3 ) ■ V^ 3 - 2Vft 3 • Vft 3 - 4ft 3 V • Vft 3 dA 

+ J fll [fc c (16iJ 4 - 20H 2 K + 4K 2 + Kcl) + 2K\i - 2HAp] dA. (26) 

The detailed expressions of operators V,V,V,V 2 ,V-V,V-Vin above equation 
can be found in Appendix D of Ref. [19]. 

Now, we discuss mechanical stability of spherical lipid membrane with radius 
R. In this case, we have H = —1/R and K — 1/R 2 . Substituting them into 
equation ©, we arrive at 

ApR 2 + 2[iR - k c c (2 - c R) = 0. (27) 

Using r = i?(sin# cos </>, sin 6* sin0, cos#) to represent the point in the spherical 
membrane, we have V = -(1/R)V, V V = -(l/i?)V 2 and V 2 = R il in8 -§g ( sin9 -§e) 
Ri ^ g -Jfrz . Under the condition (|27|l , equation l|2t)|) is transformed into 



5 2 T H = (2c fc c /i? + ApR) [ sin 6d6 [ 

Jo Jo 



d(f)Vtl 
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p7T /■27T 

+ (k c c R + 2k c + ApR 3 /2) / sin 6d9 / d0ft 3 V 2 ft 3 

Jo Jo 

Pit plTT 

+ k c R 2 / shx6d6 / d0(V 2 O 3 ) 2 . (28) 
Jo Jo 

Expand f2 3 by the spherical harmonic functions [23]: 

oo m—l 

^3 =J2 £ a* ro F; m (M), «L = (-l) m a,,_ m . (29) 

If considering V 2 Y; m = -2(Z+l)Yi m /.R 2 and f£ sin 9 d6 d(j>Y{^Y Vml = 8 mm >8 lv , 
we transform equation H28|l into 

8 2 T H = (R/2) \ai m \ 2 W + 1) - 2}{2k c /R 3 [l(l + 1) - c i?] - Ap}. (30) 
/ . m 

Obviously, 5 2 Th is a positive definite form if Ap <pi = (2k c / R 3 )[l(l + 1) — 
co-R] (2 = 2, 3, • • •). Therefore, we can take the critical pressure as 

Ap c = min{p,} =p 2 = (2fc c /i? 3 )(6 - c R). (31) 

If Ap > Ap c , the spherical bilayer will be instable and inclined to transform 
into the biconcave discoid shape. 

In this section, we have discussed the elasticity of and stability of lipid bilayer 
under the pressure. But lipid bilayer is oversimplified model of cell membrane 
which contains the membrane cytoskeleton. That is, for cell membrane, we must 
take into account the contribution of membrane cytoskeleton. 

MEMBRANE CYTOSKELETON 

In this section, we will discuss the contribution of membrane cytoskeleton to 
the free energy of cell membrane. The membrane cytoskeleton is cross-linking 
chain-like protein structure which can be thought of as a polymer membrane. 
Now we will derive its free energy by analogy with the polymer membrane [24] . 

We take de Gennes' convention [25] in the this section: the entropy S is a 
dimcnsionlcss quantity and Boltzmann factor kg is implicated in temperature T. 
If we regard the protein in membrane cytoskeleton as Gaussian chains [34], its 
root of mean square end-to-end distance is Rq ~ \/Nbo, where bo is the segment 
length of protein and N ^S> 1 is the number of segments. Assume that the 
principal radii of the membrane are much larger than i?o- If we denote the in- 
plane strain tensor by e which is assumed to be a small quantity the entropy of 
the protein chain must be the function of 2HRo, KR 2 ,, 2 J and Q because entropy 
is a dimensionless invariable under the transformation of coordinates, where H, 
K, J — tre and Q = tre are the mean curvature, the Gaussian curvature, the 
trace of strain tensor and the determinant of strain tensor, respectively. Thus 



G 



we can expand it as S ~ A X {2H R a ) + A 2 (2H R Q ) 2 + A 3 K + B 2 {2J) 2 + B 3 Q up 
to the second order terms, where A\, A 2 , A3, B 2 , B 3 are constants. There is no 
first order term of 23 in the expression of the entropy because we expect that — e 
plays the same role as e in the entropy. It is useful to write the entropy in another 
equivalent form S ~ A 2 R 2 (2H + c ) 2 + A 3 R 2 K + B 2 (2J) 2 + B 3 Q, where c is a 
constant called spontaneous curvature and we expect |coi?o| *C 1. Consequently, 
the elastic free energy per area of a membrane consisting of protein chains has 
the following form 

/ = -(Mh)TS = y[(2J) 2 - vQ] + ^(2H + c ) 2 - k'K, (32) 

where h is thickness of the membrane, M the number of protein chains per 
volume, and kd, v, k' c , k' are unknown universal constants. Here we neglect the 
entanglement of proteins. We will show kd = 4MhT and v = 1 as follows. 

Let us consider an ideal case-the planar membrane with the homogenous in- 
plane strains. In this case, H, K and cq are vanishing for planar membrane with 
symmetry between its two sides. On the one hand, equation (|32[1 is simplified 
as 

f = h [{ 2jf-vQ]. (33) 

For homogenous stain e, we can express it by its components en, £22 and e\ 2 — 
£21 = in some orthonormal coordinate system so that 2 J = en + £22 and 
Q = en£22- 

On the other hand, this structure can be compared with the structure of 
rubber. In terms of the elasticity theory of rubber [26], the deformation energy of 
a planar rubber per area can be expressed as f r — (MhT/2)[X 2 + X 2 + l/(X 2 X 2 ) — 
3], where Ai = 1 + en and A2 = 1 + £22 are extensions. For small strains, it 
is expanded to the lowest order terms as f r ~ 2MhT(e\ 1 + t\\t 22 + e 22 ) = 
2MhT[(2J) 2 - Q]. Thus we can obtain k d = AMhT and v = 1 by comparing it 
with equation 

In this section we do not discuss the mechanical property of membrane cy- 
toskeleton alone. We will discuss it with the cell membrane in the next section. 



CELL MEMBRANE WITH MEMBRANE CY- 
TOSKELETON 

The free energy of cell membrane with membrane cytoskeleton is taken as the 
sum of the free energies of lipid bilayer and membrane cytoskeleton. Thus we 
write the free energy of cell membrane under the osmotic pressure Ap as: 

T = J {k d /2)[{2J) 2 - Q]dA + J [{k c /2)(2H + c ) 2 + fx]dA + Ap J dV, (34) 

where k c = k c + k' c , cq = k c co/k c , and the term related to the Gaussian curvature 
disappears because its integration J KdA is an unimportant constant so that it 
is omitted. 
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Before taking the first order variation of the free energy (|34[) . we must in- 
troduce the strain analysis expressed by the notation of differential forms [19]. 

If a point ro in a surface undergoing a displacement u to arrive at point r, 
we have du = dr — dry and naturally Sidu = Sidr (i = 1,2,3). 

If denote dr = uj 1 ei+oj 2 e 2 and du = UiWi+U 2 w 2 with |Ui| <C 1, |U 2 | <C 1, 
we can define the in-plane strains [27]: 



£ll 



£22 



£12 



du • ei 



|dro| 
du • e 2 



Ui • ei, 



1^2=0 



|dr n 



U 2 • e 2 , 



/ du ■ e 2 
V |dro| 



du • e! 
~\dxo\~ 



1*11=0 



(35) 
(36) 



- (Ur • e 2 + U 2 • ei) .(37) 



Using Sidu = Sidr and the definitions of strains l|35[l - i|37[) . we can obtain the 
leading terms of variational relations: 



SieuUJiAuj2 = SiUJiALu 2 , 

5iE 12 LOi A ui 2 = - [wi A SiUi + Si(J 2 A lo 2 ] , 

S i e 22 uj 1 Auj 2 = ivxAdi^. 

From equations l(TB )l -l|2^ )l and ||2Hl)-(|3nil, we have [19]: 



(38) 
(39) 
(40) 



8 X T 
5 2 T 



, r ,/„ n A £ndw 2 - £i2dwi d(e 12 uj 1 + e 22 u 2 ) 
k d [-d(2J) A uj 2 - + - -]ni,(41) 

, ,,/ onA £i 2 dw 2 - £22dwi d{e lx uji + £ 12 u 2 ) 
k d {d(2J) Awi - - \n 2 , (42) 

[k c (2H + c )(2H 2 - cqH - 2K) + k c V 2 (2H) 

krl 

+A P - 2H{p + k d J) - y (as n + 2be 12 + ce 22 )]n 3 dA. (43) 
Thus the Euler-Lagrange equations corresponding to the functional (|34ll are 

kd[-d(2J) Auj 2 ~ -(s n dw 2 ~ £121^1) + -d[ex 2 u-i + e 22 uj 2 )} = 0, (44) 

k d [d{2J) Am- ^{e 12 duj 2 ~ e^dtux) - id(£nWi + £i 2 w 2 )] = 0, (45) 

Ap - 2H{n + k d J) + k c (2H + c a )(2H 2 - c H - 2K) + k c V 2 (2H) 
k d 



-(asu + 2be 12 + ce 22 ) = 0. 



(46) 



Equations l|44|) and l|45|) are called the in-plane strain equations because they 
describe the in-plane strains of cell membrane under the pressure Ap. Equation 
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146(1 is called the shape equation because it describes the equilibrium shape of 
cell membrane under the pressure Ap. 

If we take a local orthonormal coordinates (u 1 ,^ 2 ) such that the first and 
second fundamental forms of the surface are / = gn(du 1 ) 2 + 322 (dw 2 ) 2 and 
II = Lu(du 1 ) 2 + 2Li2du 1 du 2 + L 22 (du 2 ) 2 , respectively, above equations |@3J|- 
l(4U|) can be expressed as: 



hi 
kd 



(e 22 - en) 



(en - £22) 



du 2 



' du 1 
/- — d 



£22 



2e 22 ) 



'9n 



1912- 



- dei2 
du 2 
de 12 



2e v 



2ei 



Ap - 2{n + k d J)H + k c (2H + c )(2H 2 



du 2 
d\fg~22 



du 1 du 1 
c H -2K) + k c V 2 (2H) 
kd 



0, 



= 0, 



\ae n 



2&ei2 + ce 2 2] = 0. 



Obviously, if k d = 0, then equations l(4*4l and (|43)l are two identities. More- 
over equation (|4*tj|l degenerates into shape equation of closed lipid bilayers 
in this case. Generally speaking, it is difficult to find the analytical solutions to 
equations l|44 |l -l|4^ )l . But we can verify that en = e 2 2 = £ (a constant), e i2 = 
can satisfy (|44|l -146 |) for a spherical membrane with radius R if the following 
equation is valid: 



ApR 2 + (2fx + 3k d e)R + k c c (c R - 2) = 0. 



(47) 



Now we discuss the effect of membrane cytoskeleton on the stability of cell 
membrane. That is, we will calculate the second order variation of functional 
(|34l) . Additionally, we only consider the spherical membrane with uniform in- 
plane strains whose radius satisfies equation ((47(1 for simplicity. 

In Ref. [19], only the term 8\T related to the out-plane mode is calculated. 
Here we also consider the contribution of in-plane modes. Due to the notation 
of exterior differential and Hodge star *, f^i and can be expressed as £l\(x>\ + 
^2^2 = dfl + *dx by two scalar functions and x- Using equations (|15(l - (|24|l 
and ((38(1 - 140(1 . we can calculate 5fJ-, S 2 J-, S 2 J-, 5i8 2 T, 6163 J 7 , and 5 2 5 3 T from 
equations l|41 |l -l|4H |l and ((47(1 . Eventually, we arrive at 



5 2 T = 6 2 T- 



5 2 T 



5 2 T 



25 X 5 2 T + 26 1 6 3 F + 25 2 5 3 F = G x + G 2 , 



(48) 



where 



Gi = 



G 2 



n 2 3 {3k d /R 2 + (2k c c /R 3 ) + Ap/R}dA 



+ / n 3 W 2 n 3 {k c c /R + 2k c /R 2 + ApR/2}dA+ / k c (V 2 fl 3 ) 2 dA 



3fc d 
' R 



n 3 \7 2 ndA + k d J (x/ 2 n) 2 dA 

k -tj(y 2 x ) 2 dA + ^Lj x V 2 x dA. 



(49) 
(50) 
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If we take K — k c /2, K — 3k d /2, /i = fed/2, w = SI3 and \P = £1 in equations (6) 
and (7) of Zhang et aZ.'s paper [28], then G\ and G2 correspond to F\[w, ty] and 
i*2[x] m that paper under the conditions of Ap = and cqR = 2. Obviously, 
there is no coupling between modes {\} and {fi,il 3 }; but there is coupling 
between in-plane mode {il} and out-of-plane mode {^ 3 }. We will show that 
in-plane modes have quantitive effect on the stability of the cell membrane 
although they can not qualitatively modify the results of Ref. [19]. 

Because G2 is obviously positive definite, we merely need to discuss G\ . 
and Q in the expression of G\ can be expanded by spherical harmonic functions 

[23] as fl 3 = X^o£m=-i a J"»*lm(0,<A) and ^ = EtoIim=-! (, '™ } '™^^) 

with a* lm = (-l) m ai- m and b* lm = (-l) m fy _ m . It follows that 



OO I 



Gi = ^^2|a im | 2 {3fc d +[;(? + l)-2][/(? + l)fc c /i 1 ' 2 -fc c co/i?-Apii;/2]} 



;=o m=0 
00 ' 3fc 

1=0 m=0 
00 I j 

+££||[2j 2 o+i) 2 -w+i)]im 2 - ( 51 ) 

;=o m=0 

We find that if Ap < Pl = [2i(t ffi_ 1]fl + 2 MV+V^oM (/ = 2, 3, • ■ ■), then 
G\ is positive definite, i.e., the membrane is stable. We must take the minimum 
of pi to obtain the critical pressure: 



Ap c = min{pi} 



3^ + 2MyoH] < fc e [23-2 C - fl] ) (3fcdi?2 < l2lh) 

2%/ y + ^(1 - 25 E), (3fc d i? 2 > 121fe c ) 



(52) 

But if we do not consider the in-plane mode we will obtain the critical 
pressure [19]: 



^1- i 4 Y TO + ^ (2 _ CojR)) (3^2 > 16fcc) ■ ' 



Comparing equation (|52|l with (|53|) , we find that in-plane modes have quan- 
titive effect on the stability of the cell membrane although they can not quali- 
tatively modify the result without considering them. 

On the one hand, equation (|52(l includes the classical results for stability 
of elastic shells. For classic shell, the critical pressure for spherical shell is 
Ap c ~ Yh 2 /R 2 [29,30], where Y is the Young's modulus of the shell. If taking 
c = 0, kd ~ Yh, k c ~ Yh 3 and R ^> h, our equation 1(521) also gives Ap c ~ 
Yh 2 /R 2 . On the other hand, equation (|52|l includes the critical pressure for 
stability of lipid membranes. For spherical lipid membranes, the critical pressure 
is Ap c ~ k c /R 3 [11], which is the natural result of equation (|52|l only if let 
kd = in it. Otherwise, if we take the typical parameters of cell membranes 
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as k c ~ 20fc s T [6,7], k d ~ 2.4^N/m [31], ft - 4nm, R ~ 1/im, c R ~ 1, we 
obtain Ap c ~ 2 Pa from equation l|52l) . which is much larger than Ap c ~ 0.2 Pa 
without considering kd induced by membrane cytoskeleton. This result reveals 
that membrane cytoskeleton greatly enhances the mechanical stabilities of cell 
membranes, at least for spherical shape. 



In above discussions, we only consider the static elasticities of cell membranes. 
However, it is more important to understand the dynamics of cells. At least, 
we much cover the dynamic behavior of cell membrane, cytoplasm, cytoskeleton 
and nucleus. In this section, we will propose a framework of cell structure 
dynamics involving tensegrity architecture of cytoskeleton, fluid dynamics of 
cytoplasm and elasticities of cell membranes with membrane cytoskeleton. We 
make a model of a cell as follows. 

(i) Nucleus. The nucleus is thought of as a small rigid sphere in the middle 
of cell because it is a relative stiff membrane. Its inner structure is neglected. 

(ii) The tensegrity architecture of cytoskeleton [32,33]. Cytoskeleton is a 
cross-linked structure comprised of microfilaments, microtubles and intermedi- 
ate filaments. Microfilaments are better at resisting tension while microtubles 
at withstanding compression. Intermediate filaments also resist tension, but 
only for significant strains of cells. Thus the cytoskeleton is an integral system 
consists of continuum tension elements (microfilaments) and discrete compres- 
sion elements (microtubles) prestressed by osmotic pressure acting on the cell 
membranes and adhesions of other cell or extracellular matrix. The total free 
energy of cytoskeleton might be written as F cs k = F ten + F com , where F ten is 
free energy contributed by tension elements and F com comes from compression 
elements. F ten can be written as the production of tension and total length of 
tension elements for large strains because the tension is continuously transferred 
and the cntropic elasticity is not important. Its form for small strains must refer 
to the entropic elasticity of biopolymers [34] . F com depends on the stresses and 
strains in each compression elements. In fact, F cs k depends implicitly on the 
relative positions of junction points between cytoskeleton and cell membrane. 
The equilibrium conditions of force can be expressed as 



if we omit the inertial term of cytoskeleton. Here {R^} and {f;} represent the 
positions of junction points between cytoskeleton and cell membrane, and the 
forces at that junction points induced by cell membrane, respectively. 

(iii) Cytoplasm and the liquid surroundings of the cell. Cytoplasm and the 
liquid surroundings of the cell are regarded as incompressible viscous fluid. The 
dynamics might be describe by Navier-Stokes equation [35]: 



CELL STRUCTURE DYNAMICS 



fi = dF csk /dr l 



(54) 



dv/dt + (v ■ grad)v 
div v = 0, 



grad p/p + (rj/p)dxv grad v, 



(55) 
(56) 
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where v, p, p and rj are the velocity vector, density, pressure and dynamic 
viscosity of fluid, respectively. The components of viscosity stress tensor r has 
the form 

nj = rj(dvi/dxj + dvj/dxi), = 1, 2, 3), (57) 

where v\, 1)2, and V3 are, respectively, three components of velocity in three 
Cartesian coordinates x\, x%, and X3. 

(iv) The coupling between cell membrane and cytoskeleton, cytoplasm as well 
as liquid surroundings. The equilibrium equations of cell membrane under the 
interaction of cytoskeleton, cytoplasm and liquid surroundings can be expressed 
as 



kd 
2 



kd 
2 



(€22 — Cllj S 1 V <?22 TT- 



= e 1 ■ At ■ e 3 



522^j-(2eii + £22) 
y^Ji ■ eiS(r-Ti), 



'9ii 



du 2 



2e r 



dJgTi 



du 2 



(58) 



/ ,9Jgn — d 

(en — €22) n 9 v5ii"T^'( ei1 + 2622) 

= e 2 • At • e 3 + 2jfj • e 2 <5(r - r.;), 



. — dei2 d^Jlm 
'922-^— ; — ^12 



du 1 



du 1 



(59) 



Ap - 2(n + k d J)H + k c (2H + c a )(2H 2 - c Q H - 2K) + k c V 2 (2H) 
kd 

— ^-[ ae ii + Zbeu + ce 22 ] = e 3 • At ■ e 3 + • e 3 <5(r - r 4 ), (60) 



where A represent the outer quantity minus the inner one. In above three 
equations, we omit the inertial term of cell membrane because it is expected to be 
much smaller than the viscosity force. Moreover, these equations are valid only 
for small in-plane deformations of cell membranes. If we only considering lipid 
membrane, these equations are degenerated to the key equations in Ref. [36]. 

Equations (|6T]t are highly nonlinear and coupling with each other so 
that they must be solved by numerical methods. The key point is to develop a 
arithmetic to dealing with the coupling boundary conditions l|58|) - (|6U|l . 



CONCLUSION 

In this chapter, we discuss elasticities and stabilities of lipid membranes and 
cell membranes. We obtain the equations to describe equilibrium shapes and 
strains of cell membranes by osmotic pressures. We find that the critical pres- 
sure for spherical cell membrane is much larger than that of spherical lipid 
bilayer without considering membrane cytoskeleton. We also try to construct 
a framework of cell structure dynamics involving tensegrity architecture of cy- 
toskeleton, fluid dynamics of cytoplasm and elasticities of cell membranes with 
membrane cytoskeleton. It is an important direction to develop arithmetic to 
solve the coupling equations lfM|) - (|B!H) in the future. 



12 



References 

[1] Lodish, H. et al. Molecular Cell Biology; W. H. Freeman & Co.: New York, 
1999; Chapter 1. 

[2] Edidin, M. Nature Rev. Mol. Cell Bio. 2003, 4, 414-418. 

[3] Singer, S. J.; Nicolson, G. L. Science 1972, 175, 720. 

[4] Helfrich, W. Z. Naturforsch. C 1973, 28, 693-703. 

[5] de Gennes, P. G. The Physics of Liquid Crystals; Clarendon Press: Oxford, 
1975; 100-103. 

[6] Duwe, H. P.; Kaes, J.; Sackmann, E. J. Phys. Fr. 1990, 51, 945-962. 

[7] Mutz, M.; Helfrich, W. J. Phys. Fr. 1990, 51, 991-1002. 

[8] Ou-Yang, Z. C; Liu, J. X.; Xie, Y. Z. Geometric Methods in the Elastic 
Theory of Membranes in Liquid Cristal Phases; World Scientific: Singa- 
pore, 1999. 

[9] Lipowsky, R. Nature 1991, 349, 475-481. 
[10] Seifert, U. Adv. Phys. 1997, 46, 13-137. 

[11] Ou-Yang, Z. C; Helfrich, W. Phys. Rev. Lett. 1987, 59, 2486-2489. 

[12] Ou-Yang, Z. C; Helfrich, W. Phys. Rev. A 1989, 39, 5280-5288. 

[13] Ou-Yang, Z. C. Phys. Rev. A 1990, 41, 4517-4520. 

[14] Mutz, M.; Bensimon, D. Phys. Rev. A 1991, 43, 4525-4527. 

[15] Lin, Z. et al. Langmuir 1994, 10, 1008-1011. 

[16] Rudolph, A. S.; Ratna, B. R.; Kahn, B. Nature 1991, 352, 52-55. 

[17] Feng, Y. Z. Bio-Mechanics; Science Press: Beijing, 1983. 

[18] Naito, H.; Okuda, M.; Ou-Yang, Z. C. Phys. Rev. E 1993, 48, 2304-2307. 

[19] Tu, Z. C; Ou-Yang, Z. C. J. Phys. A: Math. Gen. 2004, 37, 11407C11429. 

[20] Tu, Z. C; Ou-Yang, Z. C. Phys. Rev. E 2003, 68, 61915. 

[21] An, R.; Tu, Z. C. Prerint 2003, |math-phys/0307007| 

[22] Chern, S. S.; Chen, W. H. Lectures on Differential Geometry; Peking Uni- 
versity Press: Beijing, 2001; Chapter 6. 

[23] Wang, Z. X.; Guo, D. R. Introduction to Special Function; Peking Univer- 
sity Press: Beijing, 2000. 



13 



[24] Tu, Z. C; Ge, L. Q.; Li, J. B.; Ou-Yang, Z. C. Prennt 2003, 
|cond-mat/0312319| 

[25] de Gennes, P. G. Scaling Concepts in Polymer Physics] Cornell University: 
New York, 1979. 

[26] Treloar, L. R. G. The Physics of Rubber Elasticity; Clarendon Press: Ox- 
ford, 1975. 

[27] Wu, J. K.; Wang, M. Z. Introduction to Elastic theory; Peking University 
Press: Beijing, 1981. 

[28] Zhang, Z.; Davis, H. T.; Kroll, D. M. Phys. Rev. E 1993, 48, R651-R654. 

[29] Pogorelov, A. V. Bendings of surfaces and stability of shells; Providence, 
R.I.: AMS, 1989. 

[30] Landau, L. D.; Lifshitz, E. M. Theory of Elasticity; Butterworth- 
Heinemann: Oxford, 1997; 3rd edn. 

[31] Lenormand, G. et al. Biophys. J. 2001, 81, 43-56. 

[32] Ingber, D. E. J. Cell Sci. 2003, 116, 1157-1173. 

[33] Stamenovic, D.; Ingber, D. E. Biomechan. Model Mechanobiol. 2002, 1, 
95-98. 

[34] Dio, M.; Edwards, S. F. The Theory of Polymer Dynamics; Clarendon 
Press: Oxford, 1986. 

[35] Landau, L. D.; Lifshitz, E. M. Fluid Mec/iamcs;Butterworth-Heinemann: 
Oxford, 1998; Chapter 2. 

[36] Komura, S.; Seki, K. Physica A 1993, 192, 27-46. 



14 



